home *** CD-ROM | disk | FTP | other *** search
- /*
- This file is automatically written by kwrite_info.
- */
- /* KAOS DYNAMICAL SYSTEM CLASS = class_demo */
-
- d4hammm (D4 nilpotent HAMiltonian [Minus Minus])
- A nilpotent Hamiltonian vector field which is symmetric
- under D4, the symmetry group of the square.
- The "[Minus Minus]" part of the mnemonic refers to two
- minus signs present in the expression for the vector
- field. Compare with "d4hampp."
- The vector field on R^4 = {(x,y,z,w)} is given by
- xdot = y
- ydot = x(mu-(x^2 + z^2)) + delta(x)(z^2)
- xdot = w
- xdot = w(mu-(x^2 + z^2)) + delta(z)(x^2)
- The parameters are mu and delta.
- Polar coordinates are enabled. The polar variables
- are r, rp, theta, and thetap.
-
- This system has three auxiliary functions:
- Energy = 0.5(y^2 + w^2) - 0.5 mu(x^2+z^2) +
- 0.25(x^2 + z^2)^2 - 0.5 delta(x^2)(z^2)
- AngMom = yz - xw
- t = time
-
- ----------
- d4hampp (D4 nilpotent HAMiltonian [Plus Plus])
- A nilpotent Hamiltonian vector field which is symmetric
- under D4, the symmetry group of the square.
- The "[Plus Plus]" part of the mnemonic refers to two
- plus signs which appear in the expression for the
- vecto field. Compare with "d4hammm."
- The vector field on R^4 = {(x,y,z,w)} is given by
- xdot = y
- ydot = x(mu+(x^2 + z^2)) + delta(x)(z^2)
- xdot = w
- xdot = w(mu+(x^2 + z^2)) + delta(z)(x^2)
- The parameters are mu and delta.
- Polar coordinates are enabled. The polar variables
- are r, rp, theta, and thetap.
-
- This system has three auxiliary functions:
- Energy = 0.5(y^2 + w^2) - 0.5 mu(x^2+z^2) +
- 0.25(x^2 + z^2)^2 - 0.5 delta(x^2)(z^2)
- AngMom = yz - xw
- t = time
-
- ----------
- lorenz: (LORENZ system)
- A three dimensional vector field resulting from a
- truncation of a partial differential equation which
- models two-dimensional fluid convection. The system
- was first studied by Lorenz in J. ATMOS. SCI. 63
- (1963). The vector field is defined by
- xdot = sigma (y -x)
- ydot = rho x - y - xz
- zdot = -beta z + xy
- The parameters are sigma, rho, and beta.
-
- There are three auxiliary functions:
- t = time
- x+y = sum of x and y coordinates
- x-y = difference of x and y coordinates
-
- ----------
-
-
-
-
-
- nlmathieu: (NonLinear MATHIEU equation)
- This vector field describes the motion of a
- pendulum with a periodically excited support.
- This system is described in Guckenheimer/Holmes
- pp. 29-32 and references within.
- The vector field is given by
- xdot = y
- ydot = -damp y - (omega^2 + ampl cos(t)) sin(x)
- tdot = 1
- The parameters are omega, ampl, damp.
-
- There are no auxiliary functions.
-
- ----------
- dpfosc2: (Dissipative Periodic Forced OSCillator, version 2)
- This vector field describes the motion of a damped
- and periodically forced pendulum.
- The vector field is given by
- xdot = y
- ydot = dcampl + acampl cos(2 pi f_omega t)
- - damp y - n_omega sin(2 pi x)
- The parameters are f_omega (forcing frequency),
- acampl (amplitude of "AC" forcing), dcampl (amplitude
- of "DC" forcing), damp (damping coefficient), and
- n_omega (natural frequency.
-
- This system is periodic (with period 1) in the x-coordinate.
- Time is an auxiliary function.
-
- ----------
- henonmap: (HENON MAP)
- A mapping in the plane given by
- (x, y) -> (1 + y - ax^2, bx).
- For b not zero, this mapping has an explicit inverse
- giving by
- (x, y) -> ( 1/(bx), -1 + x + (ay^2)/(b^2)
- The parameters are a and b.
-
- There are no auxiliary functions for this map.
-
- ----------
- kotorusmap: (Kim-Ostlund TORUS MAP)
- A mapping on the two-dimensional torus, ie, the unit
- square with edges identified.
- The mapping is defined by a nonlinear perturbation of
- a translation:
- f(x,y) = ( f1(x,y), f2(x,y) ) where
- f1(x,y)= x + wx - a(asymm)/(2 pi) sin(2 pi y);
- f2(x,y)= y + wy - a/(2 pi asymm) sin(2 pi x);
- An explicit inverse does not exist, but inverses may be
- computed implicitly.
- The parameters are wx, wy, a, and asymm.
-
- The variables x and y are periodic with period 1.
- This system has two auxiliary functions:
- rhox = f1(x,y) - x
- rhoy = f2(x,y) - y
- These functions are related the the concept of a
- "rotation vector." See Baesens, Guckenheimer, Kim, and
- MacKay, Preprint, 1990.
-
- ----------
- dissstandmap: (DISSipative STANDard MAP)
-
- var_label[0] = "x";
- var_label[1] = "r";
- param_label[0] = "w";
- param_label[1] = "k";
- param_label[2] = "b";
- func_label[0] = "Rhox";
- siegelmap: (SIEGEL MAP)
-
-
- param_label[0] = "rho";
- param_label[1] = "exp";
- func_label[0] = "|Rho|";
- martyd3: (MARTY golubisky's D3 symmetric map)
- This is a mapping from the complex plane to itself
- which is symmetric under the group D3. Let REAL
- and IMAG denote the operators which take the real
- and imaginary part of a complex number. Let z=x+iy
- be complex and let
- iv = alpha |z| + lambda + beta REAL(z^3).
- Then this map is defined by
- x -> iv x + gamma REAL(z^2)
- y -> iv y - gamma IMAG(z^2).
- The parameters are alpha, lambda, beta, and gamma.
-
- There are two auxiliary functions:
- Modulus = |z|
- Re(Z^3) = REAL(z^3)
-
- ----------
-
- henonheiles: (HENON-HEILES equation)
- This Hamiltonian vector field is a classic example of
- a non-integrable Hamiltonian. The system was first
- discussed in Astron. J. 69 (1964).
- We write px and py for the momenta corresponding to
- the x and y directions. Then the vector field is
- given by
- xdot = px
- pxdot = -x - epsilon 2xy
- ydot = py
- pydot = -y - epsilon (x^2 - y^2)
- The parameter is epsilon.
-
- There are three auxiliary functions:
-
- Energy = 0.5( (px)^2 + (py)^2 + x^2 + y^2) +
- epsilon (yx^2 - (1/3)y^3)
- AngMom = xy - (px)(py)
- t = time
-
- ----------
- vanderpol: (VAN DER POL equation)
- This vector field describes the motion of the Van der Pol
- oscillator. See Guckenheimer/Holmes and references within
- for a theoretical discussion of this system's dynamics.
- The vector field is given by
- xdot = y - alpha( (1/3)x^3 - x)
- ydot = -x + beta cos(omega t)
- The parameters are alpha, beta, and omega.
-
- Time is an auxiliary function.
-
- ----------
- duffing: (DUFFING's equation)
- This vector field describes a periodically forced, damped,
- nonlinear oscillator. The nonlinearity is cubic in space.
- The vector field is given by
- xdot = y
- ydot = beta x - x^3 - delta y + gamma cos(omega t)
- The parameters are delta, beta, gamma, and omega.
-
- Time is an auxiliary function.
-
- ----------
- simpletorusmap: (SIMPLE TORUS MAP)
-
- (x,y) -> ( f1(x,y), f2(x,y) ) where
- f1(x,y) = x + e(wx + cos(2 pi x) + a cos(2 pi y))
- f2(x,y) = y + e(wy + sin(2 pi x) + a sin(2 pi y))
- The parameters are wx, wy, a, and e.
-
- There are two auxiliary functions:
-
-